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Abstract 

We investigate the possibility that the modified Yang-Mills theories can produce an accelerated cosmic 
expansion. We take into account some specific non-trivial solution of the modified Yang-Mills equation 
obtained by the author earlier, which allows us to build several modifications of accelerated cosmic expansion. 
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1 Introduction 



Presently it is accepted by the scientific community that the Universe is experiencing an accelerated 
expansion. This is supported by many cosmological observations, such as SNe la [T], WMAP [2], SDSS [3] 
and X-ray |4]. These observations suggest that the Universe is dominated by the Dark Energy (DE), which 
provides the dynamical mechanism for the accelerated expansion of the Universe. Moreover, they suggest 
a* that the DE equation-of-state (EoS) parameter, 7 = — might have crossed the phantom divide 7 = — 1 from 

1 ; p 

above in the near past. In order to explain this phenomena, one can either consider theories of modified 
gravity, or field models of DE. The simplest candidate of DE is a tiny positive time-independent cosmological 
constant, for which 7 = —1. However, it is difficult to understand why the cosmological constant is about 
^f*j I 120 orders of magnitude smaller than its natural expectation (the Planck energy density). This is the so- 

called cosmological constant problem. Another puzzle of DE is the cosmological coincidence problem: why 
are we living in an epoch in which the dark energy density and the dust matter energy are comparable? 
As a possible solution to these problems various dynamical models of DE have been proposed, such as 
quintessence [5|- So far, a large class of scalar-field DE models have been studied, including tachyon [6], 
ghost condensate [7] and quintom [8], and so forth. In addition, other proposals on DE include interacting 
DE models ;9, , braneworld models [TO], and holographic DE modeles [TT], etc. The quintom scenario of DE 
is designed to understand the nature of DE with gamma across -1. The quintom models of DE differ from 
the quintessence, phantom and k-essence and so on in the determination of the cosmological evolution. 

Another class of DE models is based on the conjecture that a vector field can be the origin of DE |12].|13|. 
The YM field can be a kind of candidate for such a vector field |14].|15]. At the same time, it is well known 
that a pure YM field (with its EoS 7 = 1/3) can not provide accelerated expansion of the Universe, for 
which 7 < 1/3 is required. This is a direct consequence of conformal symmetry of the Lagrangian for a 
massless YM field. Any violations of conformal symmetry (e.g., as a result of quantum corrections [16] or of 
non- minimal coupling to gravity |17| ) give a good chance for involving YM fields in reconstruction of DE. 
The alternative method for YM fields to be involved in DE problem is consideration of some sort of modified 
YM theory |18|-|20|. In this paper, we turn our attention to the issue of the YM fields as a source of DE 
in the frame of a modified YM theory. We aimed at deriving the necessary conditions for the possibility 
of those models to explain the accelerated expansion of the Universe. We also derive the corresponding 
equations for the scale factor evolution, and briefly discuss several examples for the DE models of this kind. 

2 Equations of model dynamics 

Let us consider the following action |20j : 

S = ~ J ' d\^g[K_^ + <$>{F? k F™ k )}, (1) 

where F^ k = diW k - dkW? + f abc WiW k , and <E> is a continuously differentiable function. Variation of (UJ 
with respect to g lk yields the Einstein equation 

G lk = R lk - ^g lk R = 2 K (i 3lfc $ - 2$' F%F^ , (2) 
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where $' = d$(l)/dl, and I = F? k F aik is the invariant of the Yang-Mills fields. As it follows from the action 
{TJ, the equation of motion for the field potential turns into 



d 



(3) 



We assume that the Universe is described by a Friedmann-Robertson- Walker (FRW) geometry: 

ds 2 = N(t) 2 dt 2 - a 2 (t)(dr 2 + £, 2 (r)dQ. 2 ), (4) 

where £(r) = sin r, r, sinh r for the sign of space curvature k — +1,0, —1, consequently. To study a FRW 
solution of the equations ([2]) ,(3j) , we can directly insert this metrics into action ([T}. Whereupon we obtain 
the following effective Lagrangian density (per unit solid angle): 



At the same time, the generalized Wu-Yang ansatz for the 5*03 YM fields can be written as [21] 



ws = x 



a W(r,t) 



(5) 



W" = EuabX 



b K(r,t)-\ , /.„ x a x fl \S(r,i) 



+ S 



We can make the following substitution into this ansatz [22] : 

K(r,t) = P(r)cosa(t), S(r,t) = P(r) sin a(t), W(r,t) = d(t). 
As a result, we have the following formulae for the YM strength tensor components: 

F i = F02 = F03 = 0, F12 = e~ 1 P'(r) cos a + 1 sino 
F13 = e _1 P'(r) sin^^m sina — 1 cosa^j , F23 = e _1 sin 9^P 2 (r) — 1 



(6) 



where n = (sin 9 cos 0, sin 9 sin <f), cos 9), 1 = (cos 9 cos cj>, cos 9 sin (f>, cos 9) and m = (— sin <f), cos 0, 0) are the 
orthonormalized isoframe vectors, and the prime means a derivative with respect to r. As noted in [22], the 
YM field ([6} has only magnetic components. It is easy to find from ([4]) and (J6| that the YM field invariant 
I = F°- k F alk becomes as follows: 
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I = 



e 2 a 4 £ 2 



(7) 



Varying the effective Lagrangian density ([5]) over P(r), and taking into account (0 we obtain the following 
Euler-Lagrange equation instead of YM equation ©: 



W + p'$"— (— 

J + e 2 aUf 2 



2P 



,2 , (P 2 - 1) 



)' 



= 0, 



(8) 



where = d${l)/dl, $" = d 2 $(/)/d/ 2 . 

The nontrivial solution for the YM equation obtained in [22] . P(r) = f (r) = cos r, cosh r for = +1, — 1, 
consequently, satisfies equation ([8]). Indeed, this solution turns both additive terms in the left-hand-side of 
{8j to zero. As it follows from ([7]), the valuable feature of this solution is that the YM invariant built on it 
depends only on time: 

l = W = -£ns- (9) 



The Hamiltonian constraint for ([5]) is 



(i) 



2 k 



37rG-.,-, A 
3 v ; 3 ' 



(10) 



where 7 should be replaced with its value By variation of ([5]) over a(t) with the subsequent choice of 
the gauge N = 1, one can obtain the following Friedmann equation: 



2^ 
a 



2 + A =8 ^(/)^^<f.'(/) + A. 



(11) 



2 



3 Accelerated expansion 

First we have to note that instead of equation ()lip one frequently uses the following equation 

$(/) - 27*' (I) 



a _ 8tyG 
a 3 



which can be obtained by combining two equations (|10|) , fTT}, and with the help of ((9J. This equation 
is just the differential consequence of equation (|10p . Nevertheless, it is convenient for investigation of the 
accelerated expansion. 

Comparing equations (|10[) and (|12|l with the similar ones of the standard FRW cosmology of a perfect 
fluid, we can find the following expressions for the effective energy density and pressure of the YM field and 
the cosmological constant: 

p = HD + ^ P = -*(i) + 1/^-3^, d3) 

where the last terms are just the energy density and pressure associated with cosmological constant A with 
EoS: pa = — PA ~ „ — 77. Therefor, the EoS for the YM field and cosmological constant in our model is 

7 _ 1+ 4_T£<^ (u) 

As it follows from equation (|f 2[) or, equivalently, from the inequality 7 < —1/3 in (|14|l . the accelerating 
regime is possible if 

$(I)-2/$'(/) + -A; >0. (15) 

07TCr 

Now we are going to consider some interesting examples concerning the standard and modified YM 
theories. 
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a) For the standard YM theory $(/) = /, that is $(/) = -— . Plugging this $(/) into (fTil) 

167T 8ne- ! a 4 '(t) 

we have 



This gives the EoS 7 = 1/3 in the case of vanishing A as it must be for the pure radiation. Besides, the 
accelerating condition (|15[) turns into 

a{t) >a c = (3G/Ae 2 ) 1/4 (17) 

when A / 0. In this case, EoS (|16[) goes from 1/3 to —1 during the evolution of the scale factor a(t), and 

3 

becomes less then —1/3 as (11711 is satisfied. Plugging $(/) = into (1101) we have the following 

equation for the scale factor: 

which can be easily integrated. It should be noted that the similar equation was discussed earlier in |17| . 

b) Let us now suppose the phenomenological power-law dependence of $(/) on /: $(/) = AI n , where 
A, n are some nonzero constants. In this case, one can rewrite inequality (|15[) as 

(l-2n)A/"+^>0. (18) 

The latter means n < 1/2 in the case of vanishing cosmological constant: A = 0. At the same time, according 

4 1 
to (|14p the EoS becomes 7 = —1 + — n = constant. For n < 1/2, 7 < — — that is this model experiences 

eternal accelerated expansion. 

Let us revert to the case of non-zero A. As it follows from @ and (1141) . 

With the help of (118[1 . it is easy to show that 7„ < —1/3 in this case too. Of course, the particular case 

6™ 

n = 1, A = 1/16-7T leads to (fT6l) . Now plugging $(/) = A 2n 4n ^ m l° (HBl we have the following equation 
for the scale factor: 

•2,7 d 2(l-2n) . A 2 Cr,n\ 

a + k = Ba K ' + —a , (20) 
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where B = 8nGA6 n /3e 2n . This equation can be integrated for several n in an explicit form. 

c) Now we consider the case of a widely discussed non-Abelian Born-Infeld (BI) Lagrangian (see, e.g.. 
[17] and bibliography therein): 



W--| 



£pa paik {^P a J?o,ik^2 



•1 + 



P 2 



16/3 4 



where /3 is the critical BI field strength, F°j. is a dual YM strength tensor. From ([6]), we can find that for 
our solution the second invariant of YM field FikF' k = 0. Hence, we can identify 4>(/) with 



$(/) = 



167ra 



Vl + 2al - 1 



where a = 1/2/3 . Under our solution for YM field, this model goes from 7 = —1/3 to 7 = —1 according to 
the following equation: 



, 8a/ 4 12a\" 1 / 2 r / , 12a / A , 



(21) 



Simultaneously, the scale factor is driven by the equation 



a +k = — 
3a 



a 4 + 



12a 



a ) + —a 



(22) 



It should be noted that the same equation was obtained earlier in non-linear BI theory on the brane [17] . 

d) At last, let us consider the effective YM field cosmic model based on the effective Lagrangian up to 
1-loop order [16] . [17] . [23] . In our notation, this Lagrangian can be written as follows: 



$(/) = -^/ln 



(23) 



where k is the renormalization scale of dimension of squared mass, b = 11/127T 2 is the Callan-Symanzik 
coefficient for the generic gauge group considered here. From (|21fl , we can find that $' = -(6/4) (In |7/2k 2 | + 
1). Due to this formula together with @ and (]14[) . we obtain the following EoS: 



1 4 
3 ~ 3 



Ae 2 4 . . 
1+ 12^ a <*> 



In 



Ae 2 4 : 
+ 12^G6 a W 



(24) 



which displays more complicated behavior on time then the ones considered above. Indeed, it starts from 
7 yMC = 1/3 at a = but approaches 7 rMC = —1 through the break of its continuity at a(t) = a cr , which 
can be found from the following transcendent equation: 



Nevertheless, if the model starts its expansion at ao > a cr (non-singular model), then no problem of that 
kind occurs. Taking into account ([9]), (]10[) and (|23[) we can find that the scale factor of the model is driven 
by the following equation: 



(a) 2 + k 



4nGb 



In 



(e K ) 2 a 4 (t) 



, A 2 
+ -a . 



(25) 



In the case of vanishing A, the EoS of this model follows from (I24[l as 

(e K fa 4 (t) 



7° 
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The critical value a(t) becomes a° r = (3/e 2 k 2 ) 1 / 4 7^ 0. 

Finally, we have to note that EoS (fT9l) . (f2ll . (|24|) and corresponding equations (f20|) . (|22]| . (|25l derived 
in the sub-sections b), c) and d) consequently can be investigated in more details. That could be done 
analytically or, in any case, numerically. We do not make it our aim in this short communication. 
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4 Conclusion 



In summary, the standard and some modified YM theories in FRW cosmology are studied in this paper. 
The specific non-trivial solution of the modified YM equation JSJ proposed by the author earlier allows us 
to build several modifications of accelerated cosmic expansion. All of them possess EoS 7 ~ — 1 at late time, 
so the cosmic coincidence problem can be avoided in those models. Besides, we have derived the equations 
for the cosmic scale factor in all those models. In our opinion, more significant result of our study is that we 
can now to consider the wide range of modified YM theories in cosmology. For such a purpose, the equations 
(|9ll0p . (|13ll4p and inequality (|15[1 have to be employed. Further details and consequences of the modified 
YM models considered here are in progress. 
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